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Abstract 

We analyze, without resort to any model field-mixing scheme, the leading temperature-dependent 
term in the 'diameter' of the coexistence curve asymptotically close to the vapor-liquid critical point. 
For this purpose, we use a simple non-parametric equation of state which we develop by meeting 
several general requirements. Namely, we require that the desired equation (1) lead to correct 
asymptotic behavior for a limited number of the fluid's parameters along selected thermodynamic 
paths, (2) reveal a Van der Waals loop below the critical point, and (3) be consistent with a rigorous 
definition of the isothermal compressibility in the critical region. For the temperature interval in 
question, the proposed equation approximates experimental data with an accuracy comparable to 
those given by Schofield's parametric equation and by other authors' equations. The desired term is 
obtained by applying the Maxwell rule to the equation and can be represented as Z>2/3 |i~| 2 ^, where 
r| = \T — T c \/T c and (3 is the critical exponent for the order parameter. The amplitude D2/3 is 
determined explicitly for the volume-temperature and entropy-temperature planes. 

PACS: 05.70.Ce, 05.70.Jk, 64.10+h, 64.60.F- 
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I. INTRODUCTION 

The liquid-vapor critical point has recently attracted 
considerable interest in the context of liquid-vapor asym- 
metry in near-critical fluids [E S S 3- This asym- 
metry manifests itself as a deviation of the tempera- 
ture behavior of the "diameter" of the coexistence curve, 
d = (fli+p2)/2p c , from the critical isochore p = p c . 
Here p\ and pi are the densities of liquid and saturated 
vapor for a given temperature T, and p c is the critical 
density. The first empirical evidence of this fact was ob- 
tained by Cailletet and Mathias [5[ , who suggested a lin- 
ear behavior of d with T as the critical temperature T c 
is approached: d = 1 + Di|r|, where t = (T — T c )/T c 
is the reduced distance from T c and D\ is a constant. 
This relation is known as the "law" of the rectilinear di- 
ameter. The later theoretical investigations @, 0, and 
then [l|, 0, 0, 3 revealed certain deviations from this law, 
which can be described with two additional nonanalytic 
terms, proportional to |r| 1_Q 0,0,11 and |r| 2/3 HUSH, 
a and /? being the critical exponents for the heat capacity 
at constant volume cy and the order parameter respec- 
tively: 



d = 1 + D 



2s\r\^ + 



DAt\ 



(1) 



For fluids, 1 



> 2/3, and the first temperature- 



* Present address: Department of Advanced Mathematics, Odessa 
National Polytechnic University, 1 Shevchenko Ave., Odessa 65044, 
Ukraine. 

1 Corresponding author. E-mail address: mrs@onu.edu.ua (M. Ya. 
Sushko) 



dependent term on the right of Eq. |T]) is expected to 
dominate in the region asymptotically close to the criti- 
cal point. Since the extrapolation of d to the critical tem- 
perature is commonly used to obtain the critical density, 
the determination of this term becomes of great practical 
importance. 

Theoretically, Eq. (p} is usually justified by starting 
from scaling theory, dealing with two independent scal- 
ing fields, hi and /12, and two conjugate scaling densities, 
A\ and A2. The critical part of the thermodynamic po- 
tential is taken to be a homogeneous function of hi and 
h'2- The form of the scaling function and the values of a 
and j3 for a particular model can, in principle, be calcu- 
lated by using the renormalization group approach and 
e-expansion. For all models of the same class of univer- 
sality, these quantities are identical. In particular, fluids 
are often treated in terms of the lattice-gas model, which 
is isomorphic to the Ising model. However, this statement 
is approximate, for the Ising model has perfect symme- 
try with respect to the change of sign of the quantities 
hi and Ai, which real fluids lack. As a result, the coex- 
istence curve within this model is symmetrical (the con- 
stants D in Eq. (|TJ) are zero). In order to describe the 
critical point of real fluids and account for asymmetry of 
the liquid-vapor coexistence curve, various field-mixing 
schemes have been proposed (see Ref. [H, [l(| and further 
developments P, HJ , within which independent physical 
fields (say, the chemical potential and temperature) and 
the conjugate physical densities (the density and entropy) 
are represented as linear combinations of the scaling fields 
(hi and /12) and densities (A\ and A2) of the symmetric 
model. 

Besides solving the problems of the proper choice of 
the order parameter, the conjugate field, the form of the 
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scaling function, etc., which are difficult in themselves, 
practical implementation of field-mixing schemes implies 
the use of certain model assumptions of both the form of 
the mixing and the numerical values of the asymmetry 
coefficients in the corresponding linear combinations. In 
practice, these coefficients and, consequently, the coeffi- 
cients D in Eq. |T]) serve as adjustable parameters. The 
extraction of their numerical values from experimental 
data is a very difficult task in itself (see the discussion 
of this matter in Q), due to the difficulties associated 
with carrying out experiment in the critical region and 
the problem of separation of different terms in Eq. (Q]) . 

In view of the above-said, the question arises of 
whether there is an alternative way of justifying Eq. |T]), 
which would allow definite conclusions on the universal- 
ity of separate terms in Eq. ([T]) and estimations of their 
coefficients. In the present communication, we report the 
results of study of the leading temperature-dependent 
terms in d which were obtained without resort to any 
model field-mixing scheme. Namely, we proceed from the 
idea [ll| that an efficient equation of state for the asymp- 
totic neighborhood of the liquid-vapor critical point can 
be developed by using information on the asymptotic 
behavior of a limited number of the parameters of the 
fluid alon g se lected thermodynamic paths. The original 



approach If] is modified in several directions. Firstly, 



we expand the set of the fluid's parameters used to con- 
struct the desired equation. Secondly, we take into ac- 
count a rigorous definition of the isothermal compress- 
ibility in the critical region; this definition becomes a 
source of asymmetry of the coexistence curve. Thirdly, 
we require that for r < 0, the desired equation reveal a 
typical Van der Waals loop. Having met these require- 
ments, we obtain the desired equation and then test it 
for the accuracy of interpolation of experimental data. 
The sought-for asymptote of d is determined in the ex- 
plicit form by applying Maxwell's rule to the equation. 
Finally, we analyze the behavior of the diameter in the 
entropy-temperature plane. 



II. CONSTRUCTION OF THE EQUATION OF 
STATE 

In his work [ll[ , Martynov developed a simple equation 
of state for critical fluid by using the asymptotic law for 
the behavior of the isothermal compressibility (3t on the 
critical isochore and the asymptotic form of the equation 
for the critical isotherm. In the standard variables 7r = 
P/P c - 1, r = T/T c - I, oj = V/V c - 1 and for t > 0, 
these relations arc written as 



as oj = and r — * + , 



7r = -D oj\oj\ s 1 + ... as t = 



(2) 



(3) 



The expression for pressure obtained in Ref. [Tl| fails 
to adequately interpolate PpT data near the liquid-vapor 



critical point (see Ref. Nonetheless, its implicit 

functional form was used by the authors of Ref. [12J as 
the basis for a scaling equation of state. Invoking the 
field- mixing scheme |I0j, they developed a nonparamet- 
ric equation that interpolates experimental data with an 
accuracy comparable with that given by Schofield's para- 
metric equation (3- However, the expression for d ob- 
tained in Ref. [12[ contains only one nonanalytic term 
oc \t\ 

In the present work, we expand the original set (2), 
(3) to incorporate the well-known fact that the deriva- 
tive (dir/dr)^ remains finite along the critical isochore, 
including the critical point itself (see, for instance, Ref. 
(To|). Consequently, we additionally require that desired 
equation satisfy the relation 



da 



= M — const as r = and oj = 0. 



(4) 



Also, we pay attention to the fact that the rigorous 
definition of (3t in terms of 7r, r, and oj is written as 



'T — 77 



dV_ 



NT 



1 1 

" 1 + OJ~P c 

J_ / 01n(l-f 



du> 
dir 



(5) 



Eq. ([5]) is asymmetric with respect to the transforma- 
tion oj — > — oj, but this source of asymmetry, to our best 
knowledge, has never been given an in-depth analysis. 

With postulates we begin our consideration 

with the region t > 0. Following reasoning [Tl|, we note 
that Eq. ([3]) is actually the equation of state valid only for 
the curve r = 0. As we shift from the latter, the desired 
equation takes the form 7r(r, oj) — 7r(0, w)+/(t, oj), where 
the unknown function f(r,oj) is expected to vanish at 
t = and to satisfy Eqs. ©, flU and ©. It follows that 



7t(t, oj) = Mt 



1 



ln(l + w) r 7 -D oj\oj\ 6 



T > 0. 



(6) 

For the classical values of the exponents 7 = 1, S = 3 
and amplitudes Tq = 1/6, Dq = 3/2, and for the value 
M = 4 Eq. ^ coincides with the asymptotic form of the 
Van der Waals equation for r, oj « 1: ir = At — 6ojt — 
(3/2)oj 3 + .... 

Next, we suggest that below the critical temperature 
(t = — |t| < 0), Eq. © must be generalized so as to 
demonstrate a typical Van der Waals loop and transform, 
for the corresponding values of the exponents and ampli- 
tudes, to the asymptotic form of the Van der Waals equa- 
tion for |t|, oj « 1: 7r = -4|r| + 6w\t\ - (3/2)cj 3 + ... . 
Then 



it(t,oj) = -M\t\ + — ln(l- 



-DqojIoj 



s-i 



Combined together, Eqs. (|6|) and ([8]) give 



T < 0. 

(7) 



7t(t, oj) = Mt - 



To 



ln(l + oj)t \t\ 



7-1 



A^M^ 1 - (8) 
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Eq. {§1 represents the desired equation and forms the 
basis for further analysis. 



III. TESTING THE EQUATION OF STATE 

In the quadratic approximation in to, ln(l + uj) w 
uj — g> 2 /2, and Eq. ([8|) almost coincides with the equa- 
tion [l2|], except for terms of order I 5 " 1 " 1 and |r| 2_a . 
In its accuracy, the latter proves to be comparable with 
Schoficld's parametric equation. 

Here, we consider a more sophisticated way of testing 
Eq. ([5]) , hinted by the suggestion that the law of corre- 
sponding states can be extended to the asymptotic neigh- 
borhood of the liquid-vapor critical point. If so, then the 
variables tt, t, and ui can be properly normalized so as 
to fall on one and the same curve. For Eq. ([5]), the form 
of this curve is particularly suitable for processing, with 
only one adjustable parameter actually used. 

The presence of the logarithmic factor in Eq. ([5]) makes 
the choice of the normalized variables simple. After we 
set 



V 



TT 



t = 



(A,ro) 1/7 ' 



uj, 



Eq. |(5J) takes the form 

p(t, v) = at- In (1 + v) f - v s , 



(9) 



(10) 



where a is a constant. If the law of corresponding states 
does occur for a certain class of near-critical fluids, a 
remains the same for all members of the class. In par- 
ticular, the Van der Waals and Dieterici equations give 
estimates a = 2/3 and a = 3/4 respectively. 

To test the validity of Eq. (fTT))) for interpolation of PpT 
data in the close vicinity of the critical point, we used 
data for water 0, [Uj], nitrogen [l6[, and carbon diox- 
ide [l7| ■ The values of the critical exponents for all the 
substances were taken to be equal to their values in the 
three-dimensional Ising model: 7 = 1.239, 6 = 4.80 fl8l ]. 
The critical amplitudes were calculated by using the val- 
ues of the parameter z = P C V C /RT C and approximation 
formulas [l9[ ■ The parameters of the critical points were 
taken from the original works [HI, HE EE E3 • 

To scrutinize the functional form of Eq. (fit)]) , the val- 
ues of the function / = p + In (1 + v) t 1 — v s were 
plotted for available values of the variables ©, and 
then the possibility of their interpolation with a straight 
line at was analyzed (see Figs. 1 and 2, where Ap = 
(p — p c )/p c ). The / versus t dependence was indeed 
found to approach a linear one as the density inter- 
val was narrowed. The slope a of the interpolation 
line was determined with the least square procedure 
and then used to estimate the root-mean-square error 

/ " 2 
Af = JN ^ [(-Pcxp - -Pthcor) /-Pcxp] for direct ap- 
proximation of the PpT-data with the original equation 
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FIG. 1: / versus t in the region -2.472 x 10" 4 < r < 1.236 x 
10" 3 , -0.3 < Ap < 0.3 according to data 0, Gl| for water 
(a) and its asymptotic behavior for —0.15 < Ap < 0.15 (b); 
a = 0.63 and Af = 0.27%. 
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FIG 2: / versus t in the region -1.236 x 10" 3 < r < 3xl0~ 3 , 
—0.23 < Ap < 0.23 according to data [16| for nitrogen (a) and 
its asymptotic behavior for —0.15 < Ap < 0.15 (b); a — 0.48 
and Af = 0.4%. 



© (see Fig. 3). In all cases studied, Af was comparable 
with that obtained in Ref. 113. and even lower for water. 



IV. THE COEXISTENCE CURVE 

Let V\ and V2 be the maximum and minimum molar 
volumes (for a given r) of the coexisting liquid and vapor 
respectively, and let wi = (Vi — V c )/V c = —\^i\ and 
^2 = (V2 — V c )/V c = uj 2 > be the corresponding values 
of uj. It is well-known (see, for example, [20]) that the 
form of the coexistence curve can be determined from 
the equality of the chemical potentials for the coexisting 
phases, or, equivalently, from Maxwell's rule. 

Below the critical point, the liquid (u) < u)\) and vapor 
(uj > UJ2) phases satisfy Eq. ([8]). In particular, for the 
values ui\ and u>2 on the coexistence curve we have: 

7r(r,o;i) = -M\r\ + -^ln(l - M) |r| 7 + A> ^it\ 
l 

(11) 
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FIG. 3: Deviation (P cxp - P t hcor) /Pcx P ■ 100% of data [H] 
from Eq. ©, with a = 0.48, for isotherms r = 3.169 x 1CT 3 
(o), 2.694 x 10 -4 (+), 7.448 x 10~ 4 (x), and 3.122 x 10 -3 (o). 



tt(t,u; 2 ) = -M\t\ + + lu 2 ) |r| 7 - D lu 2 2 \uj 2 \ S 2 . 

1 o 

(12) 

For the two-phase region (uji,uj 2 ), Eq. (JSJ) must be cor- 
rected according to Maxwell's rule: 



7t(t, uj) du — tt {oj 2 — U>l), 



(13) 



where 7f is the constant pressure along the isotherm- 
isobar (lui,uj 2 ). At the endpoints of this segment, n is 
equal to the values given by (fTTjl. (fT2"j): 7f = 7r(r, a>i) = 
7r(r,wi). 

Taking the integral (|13p and combining the result with 
Eqs. (dU) and (JT2J), we find: 

^HXl + wa) -ln(l - | W1 |)] = £ + l^l 4 ) , 

(14) 

6 -D (J+I-^ 1 ' 



i |r| 7 (w 2 + |wi|) = 
1 o o 



1 



(15) 



Thus, the analysis of the coexistence curve reduces to 
the study of the system of Eqs. (fT4")l . (fTS")) . To obtain the 
shape of this curve, we define a nonnegative quantity x 
by the relation 

U2 = x\ui\, (16) 

and then divide both sides of Eqs. (fT4"]) , (fl"5)) into each 
other to eliminate the variable r. As a result, we see that 
x is the only root of the transcendental equation 

5 



1 



|wi| (x 



5+1 



1 



1 



hi 



1 + \UJl\x 



-|w 1 |(^ + l)(x + l)=0, (17) 



with |wi| being a parameter. For a fixed value of |u>i|, 
the left side of Eq. (fTT|) is a monotone increasing function 
g(x) of x. The single root of this function is close to unity 
and is located in the region x < 1. 

To obtain the analytical solution of the system {HI, 
(|15[) . we represent cc = 1+e, where e « 1. In the linear 
approximation with respect to e, Eq. (JTTJ) gives 



iKI ( 



i 



2|wi| 



hi 



l+l^l 



<f(l + M)ln± 



i+l^il 

U)l 



1+0)1 



^ (18) 

Expanding the logarithmic functions on the right of 
Eq. (|18[) in power series with respect to |o;i| and restrict- 
ing ourselves to the linear approximation, we find e = 
and oj 2 = |wi|. Eq. (fT4)) then immediately gives the well- 
known asymptotic law 



ui 2 ■■ 

with amplitude 



-wi = S |r| 



Bn 



1 



(^ r 



4/(5-1)' 



(19) 



(20) 



In the second-order approximation with respect to |ct>i|, 
these two relations still hold. Asymmetry of the coexis- 
tence curve appears only within the third-order approx- 
imation, in which case we find that asymptotically close 
to the critical point, 



c = — 



2UI 



3(5-1) 



O 



(m 2 ) , 



0J 2 - LOl = 2B \t\ 



o 



{\r\ 20 ) , 



(21) 
(22) 



lu 2 + uji — e\uo\\ 



3(5-1) 



B 2 \r\ 2 " + o(\r\ 2 ' 3 ). (23) 



It follows from Eq. (|2"3"|) that the desired asymptotic ex- 
pression for the "diameter" of the coexistence curve in 
the volume-temperature plane is 



Vt + V 2 . 1 

dv = = 1 ; 

v 2V C 3 (6 - 1 



2u|2/3 



(|r| 2/3 ) ■ (24) 



V. THE ENTROPY-TEMPERATURE PLANE 

Incorporating an additional requirement [r"fl | on the 
asymptotic behavior of the molar heat capacity cy (r, uj) 
along the critical isochore uj = in the two-phase region, 



cy(r,0) =Ai\t\ a + 



(25) 



let us determine the "diameter" d s = {s% + s 2 )/s c of the 
coexistence curve in the entropy-temperature plane. For 



5 



this purpose, we first use the relation d\x = —sdT + VdP, 
s and V being the molar entropy and volume of the sys- 
tem, to find the chemical potential /i. In terms of r, u), 
and 7r, and for a constant temperature, we have 



(dfJ,) T = V c P c (l+Uj) (dn) T . 
Correspondingly, 



H{t,u) = V" c P c tt + Lodir) +/(r), 



(26) 



(27) 



where the integral is taken along an isotherm and f(r) is 
a function of temperature alone. In view of Eq. l[8]). we 
find 



/i(r, ui) = P C V C 



Mr wsgn(r)r 2 |r 

To 



7-2 



D sgn(cj) uj^\uj\ S 2 - D " sgn (uj) uj\uj\ c 



f(r) 



(28) 



In the asymptotic vicinity of the critical point, the 
fourth term in the brackets can be neglected, and Eq. 
takes the form 



H(r,ui) - /x(r,0) 

PrVc. 



-u)\uo\ S 1 h(r I \l 



(29) 



in accordance with the scaling hypothesis [21(. The 
asymptote of the scaling function h(x) for x << 1 is 
h(x) = Do + p-.x 7 . Implicitly present in Ref. ll|, this 
form of the scaring function was later postulated in Ref. 

The omitted term in the brackets in Eq. ([28]) rep- 
resents, evidently, the first correction to the asymptotic 

law ipg]). 

As a check on the validity of Eq. (|28[) . we can easy 
verify that the phase-equilibrium conditions for two co- 
existing phases with molar volumes ui\ = — \u>i\ and uj 2 
(t < 0), 



7r(T, UJl) = 7r(T, UJ 2 ), (i(t, LUl) = (l(T, U> 2 ), 



(30) 



reduce to Eq. (|14|) and (fT5|) . provided, as was already 
suggested, that /(t) is independent of uj. 
In view of Eq. (f2"8"]) , the molar entropy, 



1 fdfi 



(31) 



s(t, uj) = ^ Smu + J- [w - (1 + w) In (1 + w)] 



c 7 sgn(r)T|r|^ 2 }-l/'(r). (32) 



Correspondingly, the molar heat capacity at constant vol- 
ume, 



is given by 



(33) 



cv(r, w) 



P C K 1 



x 7 ( 7 - l)sgn(r)(l + r) |r| 
If, according to Eq. (|25p . 

1+T 



[o;-(l+w)lii(l+a;)] 



7 " 2 -^/"W. (34) 



/"(r)=A 1 |rr« + ..., 



(35) 



then, taking into account that the critical value of the 
entropy s c = s(0, 0) = -f'(0)/T c , we find 



s(t, uj) = ^{mlo + - (1 + uj) In (1 + w)] 



X 7 



T 

sgn(r)r|rr- 2 } + Sc + | (A x |t|-« + ...) ,(36) 



whence 



T c « c 3(<J - 1) 1 1 



s c V 1 - a t c 2r 



(37) 



It follows that condition ([25)1 is not necessary to ex- 
plain the appearance of |r| 2 ^ and | t - | 1 a terms in the 



"diameter" (137)) . On the other hand, Eq. ([34)1 shows that 
on both branches of the coexistence curve, cy diverges as 
\t\ "as the critical point is approached from below. 



VI. CONCLUSION 

The temperature behavior of the "diameter" of the co- 
existence curve in the asymptotic vicinity of the vapor- 
liquid critical point has been studied within a model- 
free approach, based upon general thermodynamic def- 
initions and relationships. The critical exponent of the 
leading temperature-dependent term in the diameter is 
found to be 2(3. The critical amplitude for this term 
is determined explicitly for the volume-temperature and 
entropy-temperature planes. In the latter case, the 
"1 — a" term has been recovered as well. 
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